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1 Introduction 

In this work, we adapt approximating operators for real-valued functions to set-valued 
functions (multifunctions, SVFs), by replacing an operation between numbers with an 
operation between sets. The known approximation methods, based on Minkowski sums 
of sets, fail to approximate, when the images of a multifunction are not convex. In 
case of Bernstein-type operators and subdivision operators there is a phenomenon of 
" convexification" ( jHJ El ) . 

In [l] a binary operation between sets, the "metric average", is introduced and the 
metric piecewise linear interpolant based on it is shown to approximate continuous SVFs 
with general images. The use of this operation in the adaptation of known approximation 
methods to SVFs, requires a representation of the approximation operators by repeated 
binary averages. Such a representation exists for any operator which reproduces con- 
stants, but is not unique [TU]. This non- uniqueness leads to different operators and it is 
not clear what are the appropriate adaptations. Spline subdivision schemes represented 
by repeated averages |Hj and the Schoenberg operators defined in terms of the de Boor 
algorithm [7j are proved to approximate SVFs with general compact images. Yet, for 
the adaptation of the Bernstein operators based on the de Casteljau algorithm we could 
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obtain an approximation result only for SFVs with images in R all consisting of the 
same number of disjoint intervals 

In this paper we introduce a set-operation on a finite sequence of compact sets, 
termed "metric linear combination", which extends the metric average. We adapt ap- 
proximation methods for real-valued functions to SVFs, replacing linear combinations 
of numbers by the metric linear combinations of sets. We prove that this adaptation 
of any linear operator, approximating continuous real-valued functions, approximates 
continuous SVFs of bounded variation. In particular for Lipschitz continuous SVFs, 
sharper error estimates are obtained. Approximation results for set-valued functions 
which are only continuous, are given for a limited class of operators. It should be noted 
that our adaptation method is not restricted to positive operators. The approximation 
results are specialized to the Shoenberg spline operators and the Bernstein polynomial 
operators. Also the adaptation of polynomial interpolation to SVFs is presented as ex- 
amples of non-positive operators. This adaptation is illustrated by two metric parabolic 
interpolants. 

An outline of the paper is as follows. The next section contains basic definitions, 
notation and known results. In Section |3] we introduce the metric linear combination 
between a finite number of ordered sets and define metric linear operators for multi- 
functions based on it. In Section 0] properties of the metric piecewise linear intorpolant 
are considered. In particular a representation of it by a specific set of selections is stud- 
ied. Similar selections are used in [H] and [2] to prove the existence of a continuous 
selection for a continuous SVF of bounded variation, and of a representation of a Lip- 
schitz SVF respectively. In Section 03 we derive approximation results for the metric 
linear approximation operators, based on the results in Section |3J Finally, in Section 
we specialize these results to some classical approximation operators. 

2 Preliminaries 

First we present some definitions and notation. 

• K(R n ) is the collection of all compact nonempty subsets of R n . 

• A linear Minkowski combination of two sets A and B from K(R n ) is 



with A, fi G R . 

• The Euclidean distance from a point a G R n to a set Be K(R n ) is defined as 



\A + jiB = {\a + fib : a G A,b G 



B}, 



dist(a, B) = inf \a — b\, 



where | • | is the Euclidean norm in R n . 

• The Hausdorff distance between two sets A,B G K(R n ) is defined by 
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• The set of all projections of a £ R n into a set Be K(R n ) is 

n B (a) = {b E B : \a - b\ = dist(a, B)}. 

• For A, B e K(R n ) and < £ < 1 , the t-weighted metric average of A and B is 

A© t fi = {£a + (l-£)6: (a, b) E Tl(A, B)} (1) 

with U(A, B) = {(a,b) E Ax B : a G n A (6) or 6 G n B (a)} . 
The metric average has the metric property [1^ 

haus(A ® t B, A © s B) = | £ - s\ haus(A, B), 

haus(A© t B,A) = (1 -£)haus(A,S), (2) 
haus(A ©t 5,5) = £haus(A, J B). 

• The modulus of continuity of / : [a, b] — > X with images in a metric space (X, p) is 

u>[ a , b ](f,5) = swp{p(f(x), f(y)) : \x - y\ < 5, x, y E [a, b] }, 5 > 0. (3) 

In this paper X is either R n or K(R n ) , and p is either the Euclidean distance or the 

Hausdorff distance respectively. 

The property of the modulus that we use is 

u [aM (f,\5)<(l + \)u [aM (f,5). (4) 

• By Lip ([a, b], C) we denote the set of all Lipschitz functions / : [a, b] — > X satisfying 

p(f(x),f(y))<C\x-y\, Vx,ye[a,b], 
where £ is a constant independent of x and y. 

• A variation of / : [a, b] — > X on a partition \ = { x o < ■■■ < x^ : Xi & [a, b], i — 0, N} 
is defined by 

TV 

W,x) = ^p(/(^),/(x i _ 1 )), 

1=1 

The total variation of / on [a, b] is 

V a \f) = su P V(f,x)- 

x 

We say that / is of bounded variation if V£(f) < oo and define in this case 

v f (x) = V:(f), xe[a,b]. (5) 

It is obvious that Vf is nondecreasing. If / is also continuous then Vf is continuous as 
well. For completeness we prove it. 

Proposition 2.1. A function f : [a, b] — > X is continuous and of bounded variation on 
[a, b] if and only if Vf is a continuous function on [a,bj. 
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Proof. The sufficiency follows from 

p(f(x),f(y))<V2(f)=v f (y)-v f {x) 1 for x < y. (6) 

To prove the other direction, fix x G [a, b] and e > . By the uniform continuity of / 
on [a, b] , p(f(z), f(y)) < s/2 if \z — y\ < 5 for some 5 > . First we show that Vf is 
continuous from the left. We can always choose x — { a — x o < x i < ••• < x n = x } such 
that 

N 

K(f) < V(f, X) +e/2 = Y, P(f( x i), /(xi-i)) + e/2, 

i=i 

and x — xat_i < 5 . Thus 

N-l 

i=i 

implying that — Vf(xN_i) < e. By the monotonicity of t^we get for every 

xjv_i < y < x 

Vf{x)-v f (y) < e. 

Similarly one can show the continuity of Vf from the right. Thus we obtain that Vf is 
continuous at x and consequently it is continuous on [a, b] . □ 

From (jnj) we conclude that 

W[o,6](/,£) < U[a,b](Vf,5). (7) 

• By CBV we denote the set of all functions / : [a, b] — > X which are continuous and 
of bounded variation. 

• For a set-valued function F : [a, b] — > K(R n ) , any single-valued function / : [a, b] — ► R 
with f{x) G F(x) , Vx G [a, b] is called a selection of F . 

Definition 2.2. A set of selections of F , {f a : a G .4} , is termed a representation 

of F if 

F[x) = {f a {x) : aeA}, Vx G [a, b}. 
We denote this shortly by F = {f a : a G *4}. 

3 Linear operators on SVFs based on a metric linear 
combination of ordered sets 

In this section we introduce a new operation on a finite number of ordered sets. Using 
this operation we present a new adaptation of linear operators to multifunctions. 



4 



Definition 3.1. Let {Aq, Ai, A^} be a finite sequence of compact sets. A vector 
(a , ai, a^) with G A4 , i = 0, N , for which there exists j , < j < N such 
that 

di-x G IlA^ai), 1 < i < j and a i+1 G U Ai+1 (ai), j < i < N - 1 
is called a metric chain of {A Q , ...,A^} . 
An illustration of such a metric chain is given in Figure 13.21 




a erU o (ai) aj_i e U Aj _ 1 (a j ) cij S Aj a j+1 en Aj+1 (aj) a N e U AN (a N -x) 

Figure 3.2. 

Thus each element of each set A{ , % — 0, N generates at least one metric chain. 
We denote by CH(A , A N ) the collection of all metric chains of {A , A N } . The 
set CH(A , An) depends on the order of the sets Ai , i = 0, N . 

With this notion of metric chains we can introduce a new operation between sets. 

Definition 3.3. A metric linear combination of a sequence of sets Aq,...,An with 
coefficients X ,...,Xn G R , is 

N ( N ~\ 

\A, t = I X ^ '■ ( fl o> -> *n) e CH(A , A N ) I . (8) 

i=0 I i=0 J 

Since for two sets CH(A, B) = H(A, B) , in the special case N = 1 and A , Ai G [0, 1] , 
Ao + Ai = 1 , the metric linear combination is the metric average. The following are two 
important properties of the metric linear combination which can be easily seen from the 
definition. 

N N 

(i) XiAi = Xn-iAn., , 

N N 

(ii) For Ao, Xn such that ^ Aj = 1, XiA = A . 

i=0 i=0 

With this operation, a large class of linear operators can be adapted to SVFs. 
Let A x , x — { x 0i --j^iv} be a linear operator of the form 

N 

A x (f,x) = J2c,(x)f(x l ), (9) 
defined on real-valued functions, with domain containing x ■ 

Definition 3.4. Let F : [a, b] -> K(R n ) , x C [a, b] and let {F(xi), i = 0, N} be 
samples of F at x ■ F° r A x of the form we define a metric linear operator 

A 1 " 1 on F by 

X. 

N 

A^F(x) = A%{F,x) =®c i (x)F(x i ). (10) 

i=0 

We term this operator the metric analogue of 
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Note that due to property (ii), the metric analogue of a linear operator which pre- 
serves constants, preserves constant multifunctions. The analogue of property (ii) does 
not hold for Minkowski linear combinations with some negative coefficients, even for 
convex sets. This is one reason why only positive operators, based on Minkowski sum, 
were applied to set- valued functions. As is shown in the sequel, Definition 13.41 allows to 
define also non-positive operators. 

The analysis of the approximation properties of A^F is based on properties of the 
metric piecewise linear approximation operator. 



4 Metric piecewise linear approximations of SVFs 

^From now on F : [a, b] — > K(R n ) , { Fi — F{xi) }f =Q , where a = xq < ... < xn = b and 
X= (x ,...,iAf) denotes a partition of [a, b] . We use the notation CH = CH(F , F N ) , 
and S max = max {Si : < i < N — 1}, S m i n = min {Si : < i < N — 1} with values 
Si defined clS Si — u£^_|_i — %% ) , % — 0, ...,N — 1 . In case of a uniform partition, we 
have S max = S m i n = h = (b — a)/N and denote such a partition by xn ■ 

Definition 4.1. The metric piecewise linear approximation to F is 

S?(F,x) = {\(x)fi + (1 - \(x))f i+1 : (/ ,...,/v) G CH }, x G [a^+i], 

where 

Xi(x) = (x i+1 — x)/(x i+1 —Xi). (11) 
By construction, the set-valued function S?fF has a representation by selections 

S?F = { s( X , V ):<pe CH(F 0} F N ) }, (12) 

where s(%, <p) is a piecewise linear single-valued function interpolating the data (xi,fi), 
i = 0, N, with cp = (/ , /jv) • 

Recall the piecewise linear interpolant based on the metric average, introduced in pQ: 

S^ A (F } x) = Fi® x l(x) F i+1 , x G [ Xi , x i+1 ] 

with Xi(x) defined by (JTTJ). 

It is easy to see by the triangle inequality for the Hausdorff metric and by (j2J) that 
for a continuous set-valued function F 

haus(F(x),Sf A (F,x)) < 2uj [aM (F,S max ), x G [a, 6]. (13) 

Remark 4.2. It is not unexpected that S^ IA F = S^F . 

A. 

Indeed, for a fixed x G [xi, X{ + i] and for any y G (F, x) , 

V = K(x)fi + (1 - \i(x))f i+1 

with (f u f i+1 ) G IT(Fj, F i+1 ) . Then there exists a metric chain p = (f , / i( f i+1 , f N ) , 
if G CH , such that y = s(x,f){x) . Also it is obvious that for any x G [a, b] and any 

cpeCH, s(x,<p)(x)eS™ A (F,x). 
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In the following we show that S^F , and its piecewise linear selections 112[) " inherit" 
some continuity properties of a continuous multifunction F . The following lemma and 
corollary consider Lipschitz continuous SVFs. 

Lemma 4.3. Let F G Lip ([a, b], C) , and let x be a partition of [a,b] . Then 

S^FeL ip ([a,b},C)- 

Proof. For x, y G [a^, the claim of the lemma follows from the metric property (J2J). 

Now, let x G [xj, Xj+i] and y G [xfe, x k+ i], where 0<j<k<N — 1. Using the triangle 
inequality, (J2J) and the Lipschitz continuity of F , we get 

haus(Sf(F,x),Sf(F,y)) 

< 3 _ — haus(Fj, F j+1 ) + haus(F j+ i, F k ) H — haus(F fc , F k+1 ) 

Xj+i Xj %k+l x k 

< C(x j+1 -x + x k - x j+1 + y - x k ) < C\y - x\. 

□ 

Corollary 4.4. Under the conditions of Lemma \4.3[ for any s(x, f) in\l'2\ 

s(x,<p) e L^([a,6],£). 

The proof of this corollary is similar to the proof of the previous lemma and uses the 
observation that for k > j 

k-1 

\s(x,<p)(xj+i) ~ s(x,<p)(x k )\ < Yl \ s (x, ( p)(xi) - s(x,^)(xi +1 )\ 

l=j+l 

fc-i fc-i 

< 2^ haus(S^(F,xi),S^(F,xi +1 )) <£ 2J (xi+i - x{) = C\x k - x j+1 \. 
1=3+1 1=3+1 

Now we consider the case when F is a general continuous function. 

Lemma 4.5. Let F : [a, b] — > K(R n ) be a continuous set-valued function. Then for any 
partition x of [a, b] 

u {a>b] {S^F,5)<5u la>b] {F,6). 
Proof. By definition, for any <5 > 

u [a>b] (S?F,S) = sup{haus(Sf (F,x),S™ (F,y)) : |x - y| < 5, x,ye[a,b]}. 

In case x, y G [xj, , |x — y| < 5 , the claim of the lemma is obtained using (fTTj) . 
the metric property (J2J) and (HJ), 

I I A 

haa8(S?(F,x),S"(F,y)) = haus(F 3 -, < -^(F,^) 

/; 5 x ' (14) 



Now, let x E [xj, , y E [x k , x k +i], 0<j<k<N — 1 and \x — y\ < 5 . By the 
triangle inequality 

haus(Sf (F,x),S? (F,y)) < haus(Sf (F, x), (F, ar j+1 )) 

+ haus(5 , f(F,z i+1 ),5 , f(F,a; fc )) (15) 

while by the interpolation property of SffF and since \xk — Xj+x\ < 5 , we have 

haaB(S?(F,x j+1 ),S£(F,x k )) < u [aM (F,5). (16) 
Applying dHj) and (fTT)|) to (|15j) we obtain the claim of the lemma. □ 

Corollary 4.6. For any s(x,f) in Mty and any x,y E [xj,Xj+i] , < j < N — 1 

Wx,<P)(x) - s(x,v)(y)\ < 2u [ajb] (F, \x-y\) (17) 

Also, 

uj[ ajb ](s(x,(p),S) < Auj[ a>b ](F,8), 5 < 5 min . (18) 
The proof of this corollary is similar to the proof of assertion ()14j) . 

We cannot generalize (JTSj) for arbitrary < 5 < b — a if F is only continuous. Yet 
we can get an estimate for k>[o,6](s(x> V 9 )) ^) if ^ is continuous and of bounded variation. 

Lemma 4.7. Let F E CBV([a,b}) . Then for any s(x,<p) in FHfy . 

uJi atb] (s(x,(p),6) < 4:Uj [atb] (F,5) + U[ a fl (v F , 5) < 5u [a:b] (v F ,5). 

Proof. Denote s = s(x, (p) ■ For a given 5 > , let xE [xj, Xj+i] , yE [xk,Xk+i], 
0<j<k<N — 1, such that \x — y\ < 5 . Using the definition of s(x, f) and of S^F 
we get 

fc-i 

\s(x) - s(y)\ < \s(x) - s(x j+1 )\ + ^ \s{xi + i) - s(x/)| + |s(y) - s(x k )\ 

k-l 

<?2+i — -\ s (x j+1 ) - s(xj)\ + haus(F(a:; + i), F(x{)) + k \s(x k +i) - s(x k )\ . 

J 1=3+1 

Now, (H2J yields 

<4o; M (F,5)+y^ +i (F) < 4o; M] (F, 5) + o; M (^, 5). 
Taking the supremum over |x — y\ < <5 and using (J7J), we complete the proof. □ 
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5 Approximation by metric linear operators 



We use the metric piecewise linear approximation to obtain error estimates for metric 
linear operators. 



Remark 5.1. In contrast to our previous definition of positive operators for SVFs based 
on the metric average [?|/ ; the metric analogues hlffy of two linear operators of the 
form (0|) ; which are identical on single-valued functions, are identical on SVFs. For 
example, in J21 spline subdivision schemes are not identical to the Schoenberg spline 
operators for SVFs. 

The metric analogues of linear operators of the form (JHl), which approximate real- 
valued functions, are approximating SVFs. By (f2T)J) the approximation results 
depend on the way A x approximates piecewise linear real-valued functions. 

In what follows <fi : [a, b] x R + — > R + is a continuous real- valued function, non- 
decreasing in its second argument, satisfying <fi{x, 0) = , and S x denotes the set of 
piecewise linear continuous single- valued functions, with values in R n and knots at x ■ 

Theorem 5.2. Let A x be of the form such that for any s G S x f]Lip ([a,b], £) 



Proof. By (THl 

haus(Af (F, x), F(x)) < hans(A x (S x F, x), S™(F, x)) + haus(Sf (F, x), F(x)), (22) 
while by (JU|) 




A x (s,x) - s(x)\ < C £cj)(x, S max ). 



(21) 



Then if F 6 Lip ([a, b],C) 



haus(Af (F, x), F(x)) = 2 C5 max + CC<f)(x, 5 max ). 



haus(Af(SfF,x),Sf(F,x))< sup \A x (s{ X ,p),x) - s( X ,p)(x)\ . 



In view of Corollary 14.41 and (|2~T]) 



sup \A x (s(x,<f),x) - s(x,<p)(x)\ < C C(p(x,5 max ) 



(23) 



The proof is completed by substituting (|2~Hj) and ffTTTf) in (|2"2"j) . 



□ 
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For general continuous SVFs we cannot prove an analogous approximation result. 
Yet for continuous multifunctions of bounded variation we get a weaker approximation 
result , by applying Lemma 14.71 instead of Corollary 14.41 in the proof of Theorem 15.21 

Theorem 5.3. Let F 6 CBV([a, b}) , and let A x be of the form (0), satisfying 

\A x (s,x)-s(x)\<Cui ajb] (s,(j)(x,S max )), seS x . (24) 

Then 

haus(A* f (F,x),F(x)) = 2aj[ a>b ](F, 5 max ) + 5Cuj[ ajb ](v F ,4>(x,5 max )). 

For continuous SVFs which are not of bounded variation we can prove an approxi- 
mation result only for uniform partitions and for a limited class of linear operators. 

Theorem 5.4. Let be a linear operator of the form (0), defined on a uniform 
partition xn and let h = (b — a)/N . If 

\A N (s,x) - s(x)\ < C(j)(x,U[a,b](s,h)), seS x , (25) 

then for a continuous F 

haus(A* / (F, x), F(x)) = 2 w M (F, h) + C<f>{x, 4 u [aM (F, h)). 

The proof of this result repeats the proof of Theorem 15.21 but replaces Corollary 14.41 
by (HI) of Corollary EM 

6 Examples 

In this section we present metric analogues for SVFs of the Schoenberg spline operators 
and the Bernstein polynomial operators and give approximation results. We conclude by 
two examples demonstrating the operation of metric analogues of parabolic interpolants. 
To our knowledge so far only positive operators were applied to SVFs. The two examples 
we present assert that such interpolation between sets is reasonable. 

6.1 Metric Bernstein operators 

The Benstein operator Bjsi(f,x) for / G C[0, 1] is 

B N (f,x) = g (^(1 -*)"-*/ (j^j . (26) 

It is known (see jl], Chapter 10) that there exists a constant C independent of / such 
that 

\f(x) - B N (f,x)\ < Cu M (f, y/x(l-x)/N). (27) 
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The classical Bernstein operator for F : [0, 1] — > K(R n ) with sums of numbers replaced 
by Minkowski sums of sets is 

B%»(F,x) = £ (^jx\l - x) N ~*F (£) . (28) 

It was shown in [5| that for x G (0, 1) the limit of BjJ n (F,x) when N —>■ oo, is the 
convex hull of F(x) , therefore these operators cannot approximate SVFs with general 
images. 

In Bernstein operators for set-valued functions are defined procedurally in terms 
of the de Casteljau algorithm, with the metric average as a basic binary operation, 

Ff = F{i/N), i = 0,...,N, 
F^F l k - 1 ® 1 ^ x Ft+ 1 \ i = 0,1,..., N-k, k = l,...,N, (29) 
Bjt*(F,x) = F». 

We do not know whether these operators approximate mult if unctions with general 
compact images in R n , yet they approximate multifunctions with compact images in R 
all consisting of the same number of disjoint intervals [Zj. 

Here we investigate the metric analogue of the Bernstein operators for SVFs. 
Definition 6.1. For F : [0, 1] — > K(R n ) the metric Bernstein operator is 

B^F,x) = ^^y(l-xr-F^ 

= \lt( r l) xl ( 1 -^ N ~ t f* : (/o, ...,/*)£ raj, 

where CH = CH(F(0), F(l/N)..., F(l)) . 
By Theorem 15.21 and by (j2*7jl we conclude that 
Corollary 6.2. Let F G Lip ([0, l\,C), then 

haus(B^ (F, x), F{x)) < 2C/N + C£y/x(l - x)/N. 
Moreover by Theorem 15. HI and by (j27j) 
Corollary 6.3. Let F G CBV([0, 1]) , then 

h a us(B%(F,x),F(x)) < 2u m {F, 1/N) + 5Cu [0>1] {v F , yJx{l-x)/N). 
Since (|2l)j) does not hold for these operators, Theorem 15.41 cannot be applied. 
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6.2 Metric Schoenberg operators 



For a uniform partition xn , the "classical" set- valued analogues of the Schoenberg spline 
operators for F : [0, 1] — > K(R n ) is 



N 



S^ N (F,x) = J2^/N)b m (Nx-i) 



(30) 



i=0 



where b m (x) is the B-spline of order m (degree m — 1 ) with integer knots and support 
[0, m] , and where the linear combination is in the Minkowski sense. An example, given 
in jHj, shows that the operators in (|30|) with m = 2 and iV — > oo cannot approximate 
F with general compact images, in any point of [0, 1] \ xn ■ 

A Shoenberg operator based on the metric average is introduced in j2j , by a procedu- 
ral definition in terms of repeated binary averages according to the de Boor algorithm. 
It is proved that for Holder continuous set-valued functions, the approximation rate is 
the Holder exponent. 

Here we consider the metric analogue of the Schoenberg operators. 

Definition 6.4. The metric Shoenberg operator of order m for a set-valued function 
F : [0, 1] — > K(R n ) and a uniform partition xn is defined by 



A' 



i=0 



N 



S% N (F, x ) = Qb m (Nx-i)F[-)=\J2 b ™( Nx - i )f*-- (/o. - M G CH 



i=0 



where CH = CH(F(0), F(l/N)..., F(l)) . 

By Theorem 15.41 and the known approximation result in case of single- valued func- 
tions (see [3], Chapter XII), we obtain 



Corollary 6.5. Let F be a continuous SFV defined on [0, 1] . Then 



haus(5™ JV (F, x), F(x)) = 21 + 2 



m + 1 



u m {F,l/N), xe 



m — 1 



with \t\ the maximal integer not greater than t . 



The approximation result in the specific case of Lipschitz continuous SVFs, can be 
further improved by applying Theorem 15.21 



Corollary 6.6. For F G Lip ([0, 1], C) , 

haus(^(F,x),F(x)) 



2 + 



m + 1 



C 
N' 
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6.3 Metric polynomial interpolants 

Definition 6.7. Let F : [a, b] — > K(R n ) , and let \ be a partition of [a, b] . The metric 
polynomial interpolation operator is given by 

N ( N ~\ 

P?{F,x) = ®li{x)F{xi) = <Y,li{x)fi ■ (fo,..,f N ) E CH(F(x ),...,F(x N )) \ , 

i=0 I i=0 ) 

with h(x) the i -th Lagrange polynomial, 

N 

h( X )= n 

A. A. <y . ry> . 

To illustrate our method we apply the metric parabolic interpolation operator to 
three sets in R . We consider two different examples. 
The first example: x = , X\ = 2 , x 2 = 6 ; 

F(x ) = [2, 8], F( Xl ) = {5}, F(x 2 ) = {5}. 

The second example: x = , x\ = 4 , x 2 = 8 ; 

F(x Q ) = [2,4] U [6,8], = [4.5,5.5], F{x 2 ) = [2, 4] U [6, 8]. 

The two set- valued interpolants are illustrated in Figure and Figure I6~TH respectively. 




Figure 6.9. Parabolic interpolation - second example. 

In the above figures the sets in black are F(xq) , F(xi) , F(x 2 ) and the gray curves are 
the parabolic interpolants to the selections in (|12jl. 
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